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Abstract 

Geometric properties of operators of quantum Dirac constraints and physical 
observables are studied in semiclassical theory of generic constrained systems. The 
invariance transformations of the classical theory - contact canonical transforma- 
tions and arbitrary changes of constraint basis - are promoted to the quantum 
domain as unitary equivalence transformations. Geometry of the quantum reduc- 
tion of the Dirac formalism to the physical sector of the theory is presented in the 
coordinate gauges and extended to unitary momentum-dependent gauges of a gen- 
eral type. The operators of physical observables are constructed satisfying one-loop 
quantum gauge invariance and Hermiticity with respect to a physical inner prod- 
uct. Abelianization procedure on Lagrangian constraint surfaces of phase space is 
discussed in the framework of the semiclassical expansion. 

1. Introduction 

In this paper we discuss geometric properties of the quantum dynamical systems sub- 
ject to first class constraints [[I], |^, |3|]. An importance of this problem arises in view of 
the growing interest in quantization of the so-called zero modes or moduli in various con- 
texts of high energy physics ranging from theory of gauge fields to quantum cosmology 
and nonperturbative D-brane aspects of superstrings. In contrast with oscillatory de- 
grees of freedom, quantization of these modes is impossible in a conventional Fock space 
but, rather, demands the coordinate representation of canonical commutation relations 
[|J. Quantum operators represented on configuration space of the theory as differential 
operators have interesting geometric properties with respect to diffeomorphisms of con- 
figuration space and group manifolds and, thus, deserve careful analysis for the purpose 
of consistency, uniqueness of quantization, and other subtle issues. 

At the classical level systems subject to first class constraints are described by the 
canonical action of generic form || 

S = J dt{p i q i -H (q,p)-\»T fl (q,p)} (1.1) 
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in the configuration space of coordinates and momenta (q,p) = {q l ,Pi) and Lagrange 
multiplyers A M . The variation of the latter leads to the set of nondynamical equations - 
constraints 

T fl (q,p)=0. (1.2) 

The first class constraint functions T^(q,p) and the Hamiltonian H (q,p) satisfy the 
Poisson-bracket algebra 

{T II ,T V } = U£,T X , (1.3) 
{H Q ,T U } = U£ V T X , (1.4) 

with the structure functions = U^ u (q,p) and Uq v = UQ V (q,p) which can generally 
depend on phase-space variables of the theory. 

The first class constraints indicate that the theory possesses a local gauge invariance 
generated in the sector of phase-space variables by constraints themselves T^(q,p) and 
by certain transformations of Lagrange multiplyers || . The dimensionality of the gauge 
group coincides with the dimensionality of the space of constraints, the both being enu- 
merated by the gauge index /i. If we denote the range of index i by n, i — 1, ...n, and that 
of jj, by m, \i — 1, ...m, then the number of the physical dynamically independent degrees 
of freedom equals n — m: 2(n — m) physical phase-space variables originate from the initial 
In variables {q l ,Pi) by restricting to (2n — m)-dimensional constraint surface (|1.2|) and 
then factoring on this surface out the action of m gauge transformations generated by T M . 



Dirac quantization of the theory (|1.1 ) consists in promoting initial phase-space vari- 
ables and constraint functions to the operator level (q,p,T^) — > (q,p,T^) and selecting 
the physical states | V^) in the representation space of (q,p, T M ) by the equation 

T M |*) = 0. (1.5) 

Operators (q,p) satisfy canonical commutation relations [q k ,pi] = ihSi and the quantum 
constraints as operator functions of (q,p) should satisfy the correspondence principle 
with classical c-number constraints and be subject to the commutator algebra 

[f ll ,f v ] = ihU^f x . (1.6) 

with some operator structure functions standing to the left of operator constraints. 
This algebra generalizes ( |1.3| ) to the quantum level and serves as integrability conditions 
for equations (|1.5|). 

Classically the theory ( |1 . 1| ) and this reduction to its physical sector has two types of 
invariances: the invariance with respect to canonical transformations of the initial phase- 
space variables and the geometric invariance with respect to the transformations of the 
basis of constraints. The latter property means that one and the same constraint surface 
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(|1.2|) is determined not just by one specific choice of the set of constraint functions T^q, p) , 
but by the equivalence class of those differing from one another by linear recombinations 



T; = fi^T„ n; = n;(q, P ), detn^o (1.7) 

with arbitrary invertible matrix function of (q,p) acting in the vector space of gauge 
indices. A natural question arises whether these invariances can be preserved also in the 
Dirac quantization procedure? 

The program of finding such quantum constraints for a generic constrained system has 
been partly implemented in || in the lowest nontrivial order of semiclassical expansion 
in h. The symbols of operators (T M , U^ u ) with linear in h quantum corrections have been 
found and a partial answer was given to the question of the above type: the obtained 
operators turned out to be covariant with respect to contact canonical transformations of 
the initial phase space, provided the Dirac wavefunctions ( q \ *Sf) — *f?(q) in the coordinate 
representation of commutation relations for (q,p) transform as 1/2- weight densities on the 
configuration-space manifold of q l . 

It turned out, however, that the operator algorithms for (T M , U^ u ) involve not only 
their classical counterparts featuring in (|1.3| ) but also the higher-order structure functions 
of the canonical gauge algebra @, 0. Generally the gauge algebra involves the whole 
hierarchy of structure functions and relations which begin with T^(g,p) and (|1.3|) 

G={T^ V ,U^ X ,-} (1-8) 

and at any new stage iteratively build up as consistency conditions for those of the previous 
stages. For example, the cyclic Jacobi identity {T^, {T a , T x }} + cycle(/i, a, A) = applied 



to [173] results in the equation 

{T M , U: x } + U^Ufr + cycled, a, A) = Uf uX T Pl (1.9) 

multiplied by T a , necessarily generating a new structure function U^ x antisymmetric in 
upper (and lower) indices 0. For constraints forming the closed Lie algebra all higher- 
order structure functions are vanishing, but this property depends on the choice of basis 
of constraints: the rotation of the constraint basis (|1 . 7|) can convert the Lie algebra (even 
Abelian one with U" x = 0) into an open algebra with the infinite set of structure functions. 
Thus the invariance of the the theory with respect to transformations of the form ( |1 . 7| ) 
with arbitrary Q 1 ^ necessitates considering higher-order structure functions ( |1.8| ) and their 
operator realization G — > G = {T^, C/° , U^ x , ■■■}■ 

In this paper we show that the operators constructed in really possess the expected 
properties of invariance with respect to the transformation of the constraint basis ( |1 . 7|) 
- the result briefly reported earlier in ||. Since we restrict ourselves with the one-loop 
(linear in h) approximation, we shall focuse at the covariance of only the Dirac constraints 
T M and Dirac equations on physical states ( p..5|) : the higher-order structure functions 
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will not be important for us because they are responsible for multi-loop orders of the 
semiclassical expansion. 

The paper is organized as follows. In Sect. 2 we extend the operator realization of 
quantum constraints of || to physical observables. Sect. 3 contains the proof of the co- 
variance of equations (|1.5|) , which induces the weight properties of Dirac wavefunctions in 
the space of gauge indices. These properties guarantee the unitary equivalence of quan- 
tum theories starting with different choices of constraint bases, equipped with a correct 
physical inner product. Sect. 4 gives the the fundamental two-point solution of quan- 
tum Dirac constraints explicitly featuring the weight properties of the above type. The 
reduction of this solution to physical sector is performed in Sect. 5, where this reduction 
(previously known only for a narrow class of coordinate gauges |13|, ||, |J) is extended to 
unitary gauges of general type. In Sect. 6 we consider a number of other issues omitted 
in paper |J: operator gauge independence of matrix elements of physical observables and 
their Hermitian conjugation properties in the physical inner product. The concluding 
section contains a brief discussion of possible applications of these results, while two ap- 
pendices give a number of technical details regarding the Hermiticity of observables and an 
abelianization procedure for semiclassical constrained systems on Lagrangian manifolds 
of phase space. 

2. Operator realization of quantum constraints and 
physical observables 

The operator realization of quantum Dirac constraints and lowest order structure 
functions was found in || in the form of the normal gp-ordering of their gp-symbols 
expanded up to the linear order in h. This representation implies that for any operator 
G = {Tp, U^ u , Uwx, •••} one can put into correspondence its normal symbol - a c-number 
function on phase space G(q,p) - such that the operator G can be obtained from G(q,p) 
by replacing its arguments with noncommuting operators with all the momenta standing 
to the right of coordinates. For a symbol expandable in momentum series 

oo 

G(q,p) = J2G il - in (Q)Ph-Pin (2-1) 

n=0 

this means that 

oo 

G = K P G(q,p) = E G h - in (Q)Pn-Pi n - (2-2) 

n=0 

The one-loop (linear in h) algorithms of 0] for two lowest-order operators G have the 
form 

^ = M r '-TS + T^ + < a1 )}- (2 ' 3> 
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u, 



ihd 2 U, 



fJ,V 



2 dq^dpi 



[2A) 



involving, as it was mentioned in Introduction, the higher-order classical structure func- 
tions Up£ a . As shown in [[J, these operators have two important properties. First, they 
are covariant under contact canonical transformations of (q, p) 



q l = q\q') } p, b = p k 



dq 



k' 



dq 1 



G(q,p) = G'(q',p' 



(2.5) 



under which the constraints and structure functions ( |1.8|) (all quantities bearing only 
gauge indices) behave like scalars. In the coordinate representation of canonical commu- 
tation relations the covariance of operators fl2.3| )- fl2T4|) can be written down as 



dq' 



dq 



-1/2 



G 



dq 1 



dq 



1/2 



G\ 



(2.6) 



where the operators G' are constructed by the above algorithms from their primed classical 
counterparts. This transformation law obviously implies that the Dirac wavefunction 
satisfying quantum constraints (|1.5|) should be regarded a scalar density of 1/2-weight 



*(<Z) 



dq' 



dq 



1/2 



(2.7) 



in complete correspondence with the diffeomorphism invariance of the auxiliary inner 
product of unphysical states 



{21 



This inner product diverges for physical states satisfying quantum Dirac constraints and, 
therefore, plays only an auxiliary role. It appears as a truncated inner product in the 
bosonic sector of the extended configuration space of the BFV (BRST) quantization [|7|, 
^, |9], [TO], |TT| and also will be used below (in Sect. 6) for the construction of the physical 
inner product by means of a special operator measure. 

The second important property of these operators with respect to the inner product 
(|2.8|) is their anti-Hermitian part. It is given by the trace of the structure functions and 
for Dirac constraints has the form 

%-fl = th(U^ + 0(h 2 ). (2.9) 

The algorithms ( |2.3| )- (|2T4] ) were derived in H solely as a solution of the commutator 
algebra ( O ). This helps to extend these algorithms for obtaining another class of op- 
erators - the operators of physical observables. Classically the physical observables Oi 
(enumerated by some index /) are defined as a functions on phase space, invariant under 
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the action of canonical gauge algebra. This invariance generally holds in a weak sense, 
that is only on the constraint surface 



{Oj,T,} = U^T x , (2.10) 

the gauge transformation of Oj being a linear combination of constraints with some 
coefficients Uj = Uj (q,p). (H (q,p) and Uq with eq.(|1.4j) present the example of 
such an observable and its weak invariance.) Note that again due to the rotation of the 
constraint basis (|1.7|) we have to consider nonvanishing coefficients Uf which can always 
be generated even for strongly invariant observables by a transition to another basis of 
constraints. 

In addition to their weak invariance ( 2.10|) we shall assume that the classical observ- 



ables commute with one another or form a closed Lie algebra in a weak sense 

{O h Oj} = U^jO L + U}jT x , Ufj = const. (2.11) 

In this case, from the viewpoint of commutator algebra the physical observables do not 
differ from constraints. The only difference is that unlike constraints they are not con- 
strained to vanish. Therefore, to promote the classical observables to the quantum level, 
(Oj, Uj ) — > (Oj, Uj ), and enforce the quantum gauge invariance of their operators 

[d J ,f (t ]=ihU} li f x , (2.12) 

one can use the algorithm analogous to (|2.3|) solving this commutator algebra 

6 ' - »- {°< - f wk + + j 11 " + 0{R2) } ■ (2 - 13) 

with higher-order structure functions Uj of the classical algebra (p . 1 0|) and ( |2.11|) (deriv- 
able by the method mentioned in Introduction). 

The quantum observables (|2.13|) solve the closed commutator algebra (|2.12|) . The proof 



of this statement goes by collecting the observables together with constraints into one set 
and repeating the derivation of || . The only thing to check is if the resulting commutator 
algebra does not contain nonvanishing components Uf of the operator structure functions 
(violating the weak quantum gauge invariance of observables ( 2.12Q ). This component can 



get a nonvanishing contribution only due to a higher-order structure functions Uf^ x + UfJ^ 
of the classical algebra ( ^.10[ ) and ( [2.111) . But it is easy to show that for a closed Lie algebra 



(12.11 ) the nonvanishing components of the second-order structure functions cannot have 
nongauge upper indices {U JX = 0, U JL = 0) and, therefore the quantum observables 



remain weakly invariant. For the same reason eq.( 2.14j ) does not involve the contraction 

.j 

n J - 



U? J 
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The quantum observables like constraints have anti-Hermitian part with respect to the 
auxiliary inner product ( |2.8| ). It is given by two contractions of structure functions Uj X 
and Ujj. For all reasonable compact groups generating algebras of observables the latter 
is vanishing Ujj = 0, but Uj~ x is generally nontrivial, depends on the choice of constraint 
basis and violates Hermeticity of observables in the auxiliary inner product. It is however 
inessential, because only the physical inner product must generate real expectation values 
of observables, and this will be shown to be true below. 



3. Quantum transformation of the constraint basis 

Under the linear transformation of the classical constraint basis ( |1.7|) the structure 
functions transform as 



so that 

U% = ^U x aX + {^,TJ - {lnO,^T a }, 
Q = det 02. 



(3.1) 

(3.2) 
(3.3) 



The quantum constraints T' based on the transformed basis of classical constraints 
(|1.7|) and structure functions (|3.1j ) take on the use of the algorithm (|2.3|) (with primed 
quantities) the form 



7" = Af 



^„-«^|f-f{lnii,^f„} + 0(R 2 ) 



(3.4) 



where T„ is a normal gp-symbol of constraints in the original basis and 



M " 2 dq k dp k 



^ + o(h 2 ), n = detn 



(3.5) 



From the two simple identities for operators and their normal gp-symbols 

dF l dF 2 



F 1 F 2 =M 



FiF 2 - ih 



dpk dq h 
ih 



+ 0(h 2 



F--{]nQ,F} + 0(h 2 



it is then easy to find the final form of the transformation law for quantum constraints 
under the transformation of their classical basis (17?) 



Q, = det fifj. 



(3.6) 
(3.7) 
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Similarly to ( |1.7|) this transformation involves a linear recombination of constraint 
operators with operator- valued matrix Cl" (standing to the left of constraints). This 
matrix is obtained from its classical counterpart Q^(q,p) by the algorithm Q3.7|) similar 



to ( |2.3| ): its symbol involves analogous quantum corrections except the anti-Hermitian 



part. In addition to linear combinations the transformation (3J3) includes the canonical 
transformation generated by the square root of its determinant (l 1 ^ 2 . This canonical 
transformation implies that the physical states satisfying Dirac constraints ( |1.5|) transform 
contragrediently to (|3] 



l*)'^- 1 / 2 !*) (3.8) 
and turn out to be scalar densities of weight —1/2 in the space of gauge indices. This 



property has been observed for systems subject to constraints linear in momenta in |12 
and, as we see, turns out to be true for a generic case at least in the one-loop order of the 
semiclassical expansion. 

It is easy to repeat now similar calculations for the operators of physical observables 
(|2.13|). In view of the transformation law for the trace of the observable structure function 



Uj\ = Uj X — {lndetfi, 0{\ these calculations immeadiately show that the operators of 
observables also transform canonically 

6'j = h-^Ojh 1 ' 2 . (3.9) 

Obviously, the theories differing by the choice of constraint basis should be unitarily 
equivalent. This means that the physical inner product of states | should contain 
a measure depending on this choice and transforming contragrediently to ( |3.8| ). In the 
next section we show that semiclassical states and their inner product really satisfy the 
properties compatible with the transformation law of the above type. 



4. Semiclassical physical states 



Semiclassical expansion of the operator symbols of the above type makes sense when 
the corresponding quantum states also have a semiclassical form. In the coordinate rep- 
resentation semiclassical wavefunctions 



= P(q) exp 



h 



S(q) 



(4.1) 



are characterized by the Hamilton- Jacobi function S(q) and preexponential factor P{q) 
expandable in ^-series beginnig with the one-loop ordeiQ 0(h°). The action of the operator 

lr rhis corresponds to the fact that the tree-level part is entirely contained in the exponential and is 



O^" 1 ). When the Hamiltonian Ho in eq.(Ll) is nonvanishing the wavefunction (4J), its Hamilton- 
Jacobi function, two-point kernel, etc. are time-dependent. In what follows we shall, however, omit the 
time label, because we will be mainly interested in constraint properties rather than the dynamical ones. 
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F on such functions reads as 
F*(g) 



-Pexp 



ft 



(4.2) 



where, as in (|2.2|) , F is a normal gp-symbol of F. 

The general semiclassical solution of quantum constraints (|1.5| ) with operators ( p.3|) 
was found in |13, |l| [lj]] in the form of the two-point kernel K(q,q') "propagating" the 
initial data from the Cauchy surface throughout the whole superspace of q 



K(q, q) = P(q, q) exp 



-^S(q,q') 



(4.3) 



In both expressions 
Jacobi equation 

as 

dq 



T 

-L II 



11) and (|4.3|) the phase in the exponential satisfies the Hamilton- 



(4.4) 



while the one-loop preexponential factor is subject to continuity type equation originating 
from the full quantum constraint in the approximation linear in fi 



8 

— (V* P 2 

dT 



K x p\ 



p = dS /dq 



(4.5) 
(4.6) 



For a two-point kernel the Hamilton- Jacobi function coincides with the principal 
Hamilton function S(q, q') (action on the extremal joining points q and q') and the solution 
of the continuity equation can be found as a generalization of the Pauli-Van Vleck-Morette 
ansatz for the one-loop preexponential factor [15|, 0] of the Schrodinger propagator. This 
generalization is nothing but a Faddeev-Popov gauge-fixing procedure for a matrix 
of mixed second-order derivatives of the principal Hamilton function 



>ik' 



d 2 S(q,g') 
dq 1 dq k ' 



which is degenerate in virtue of the Hamilton- Jacobi equations 
zero-value eigenvectors (|4.6j) and analogous right zero-vectors [IT 



vis 



ik' 



0, s ik ,v k , = 0, V 



dT u (q>,p') 



dp' k 



P 



-dS/dq' 



(4.7) 



giving rise to the left 



(4.8) 



Another way to view this is to parametrize time and the conjugated Hamiltonian Ho(q,p) among the 
phase-space variables and regard the Hamilton-Jacobi and Schrodinger equations as one extra classical 
and quantum constraint correspondingly [^[ [li] . 
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The preexponential factor reads 
detF ik , 



J(q)J(q') detc 



1/2 

(4.9) 



where is a nondegenerate matrix of the initial action Hessian ( |4.7|) supplied with a 
gauge-breaking term 

F w = S ih , + x?c ia ,)&, (4.10) 
and J(q) and J(q') are the Feynman-DeWitt-Faddeev- Popov "ghost" determinants JT7, 



r^] compensating for the inclusion of this term. The gauge-breaking term and ghost 
determinants are constructed with the aid of two sets of arbitrary covectors (XiiXk') 
("gauge" conditions) satisfying the only requirement of the nondegeneracy of their ghost 
operators [13|, [14], |j 



JS{q) = Xt Vt, = detJ^(g) ^ 0, 

■W) = X?Vj, J(g') = detJ^(g') ^ 0. (4.11) 

The invertible gauge- fixing matrix and its determinant (contribution of Nielsen-Kallosh 
ghosts) are the last ingredients of the generalized Pauli-Van Vleck-Morette ansatz ( |4.9| ). 

Notice now that under the transformation of the basis ( |1.7| ) the vectors ( |4.6| ) defined 
on the Lagrangian manifold of phase-space p = dS/dq transform covariantly with respect 
to their gauge indices (V^)' = ^V^. Therefore the ghost determinants transform as 
densities 

J' = (det fi£) J, (4.12) 

whence it follows that the two-point kernel with respect to each of its arguments trans- 
forms in accordance with the law ( |3.8|) in which the action of the operator f2~ 1//2 semi- 
classically boils down to the multiplication with [detf^(g, dS/dq)] . 



5. Reduction to the physical sector: unitary gauge 
conditions of a general type 

The interpretation of the semiclassical state ( (4.3|) , ( |4.9| ) is rather transparent in the 
physical sector of the theory [|T^, |5|, [TJJ. The sector explicitly arises after the reduction 
to physical variables by disentangling them from the original phase space of (q\Pi) in 
a unitary gauge. Such a reduction in |13|, ||, [14| was given for a special type of gauge 



conditions imposed only on phase space coordinates q l . Here we generalize this reduction 
procedure to unitary gauges of a general type. 
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The theory with the action ( p. . 1|) is invariant under the set of canonical transformations 
of phase space variables generated by constraints and special (noncanonical) transforma- 
tions of Lagrange multiplyers |3|, [L4| 

F<t = {<t,T ll }r», ^ = {9',T,}F, (5.1) 
5 r x n = jrn_ U^ X \ V T X - Ug x F x (5.2) 

with arbitrary infinitesimal gauge parameter JF^ = T^{t). This invariance implies that the 
equivalence class of variables belonging to the orbit of these transformations corresponds 
to one and the same physical state. The description of this state in terms of physical 
variables rather than the equivalence classes of original degrees of freedom consists in 
singling out the unique representative of each such class and treating the independent 
labels of this representative as physical variables. This can be attained by imposing on 
original phase space variables the gauge conditions 

X ^(q,p) = (5.3) 

which determine in the 2n-dimensional phase space the (2n — m)-dimensional surface 
(remember that n is the range of index i, while m is that of fi) having a unique intersection 



with the orbit of gauge transformations (5.1). At least locally, the latter condition means 



the invertibility of the Faddeev-Popov matrix with the nonvanishing determinant 

J(q,p) = detJ^q,p), J£(q,p) = { X ",T V }. (5.4) 



Gauge conditions of the form ( |5.3| ) are called unitary. They impose restrictions only on 
phase space variables and, thus, locally in time single out the physical sector and allow 
one to formulate the gauge theory in a manifestly unitary form. In this respect they differ 
from nonunitary gauge conditions imposed on Lagrange multiplyers A'' (together with 
(q,p)) and/or their time derivatives R ITJ . 



5.1. Coordinate gauge conditions 

Unitary gauge conditions take the simplest form when they are imposed only on phase 
space coordinates, x M (<?) = 0. Such coordinate gauge conditions determine the embedding 
of the (n — m)-dimensional space £ of physical coordinates directly into the space of 
original coordinates q % - superspace. This fact strongly simplifies the reduction of the 
semiclassical kernel (|4.3| ), ( |4.9| ) to the physical sector, because this reduction in the main 
boils down to the embedding of the arguments of K (q, q') into the physical subspace E. 
The geometry of this embedding, considered in much detail in JT4" |, can be better described 



in special coordinates on superspace q 1 = (£ A , 9^), in which £ A , A = 1, ...n — m, serve as 
intrinsic coordinates on £ (physical configuration coordinates), and 9^ is determined by 
gauge conditions: 

q l - q 4 = (£ A , n <f = ^(e 4 , n = X"{q)- (5.5) 
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The equation of the surface E in the new coordinates is 6^ = 0, so that its embedding 
equations coincide with the above reparametrization equations at 9^ = 0, e l (£) = e*(£, 0) 

T.:q i = e%), X^iO) = 0- (5.6) 
The relation between the integration measures on superspace dq = d n q and on E, dl; = 

d£ = dq6( X )M, M = (det[4e;])-\ (5.7) 

involves the Jacobian of this reparametrization, built of the basis of vectors tangential 
and normal to E: 

e i A = de i /d£ A , e ; = <9e7<9r. (5.8) 
Note that m covectors normal to the surface can be chosen as gradients of gauge conditions 

Xf = ^, tf4 = #, (5-9) 

that can be identified with auxiliary covectors participating in the algorithm for the preex- 
ponential factor (|4.9|) . With this identification the Faddeev- Popov operator Jtf(q, dS /dq) 
coincides with the operator J^{q) of this algorithm (which explains the use of the same 
notation). 

On the same footing with (e A , e l ) as a full local basis one can also choose the set 



(e A , V* ) with vectors transversal to E given by eq. (O ). The normal vectors of the 
first basis when expanded in the new basis 

4 = J~l V Vi + Ufa (5.10) 

have one expansion coefficient always determined by the inverse of the Faddeev-Popov 
matrix J~ x " and, thus, independent of the particular parametrization of E by internal 
coordinates. The second coefficient is less universal and depends on a particular choice 
of this parametrization. Missing information about Q A does not prevent, however, from 
finding the relation between the determinants of matrices of the old and new bases 

det [e* A , V;] = ^. (5.11) 

The reduction to physical sector in coordinate gauges follows after identifying £ with 
the physical coordinates. The corresponding conjugated momenta ha can be found from 
the transformation of the symplectic form restricted to the physical subspace (|5.6|) 

dt v 4 = Jdt [p ie \i A + p^^ : %^ ] j , (5-12) 

iiA=Pie A , (5.13) 



12 



as projections of the original momentum to the tangential components of the basis (5. 
The normal projections of Pi should be found from constraints ( |1.2| ), the local uniqueness 
of their solution being granted by the nondegeneracy of the Faddeev-Popov determinant. 
Together with ( |5.6| ) this solution yields all the original phase space variables (q l ,Pi) as 



known functions of the physical degrees of freedom (£ a ,7Ta)- The original action ( |1 . 1| ) 
reduced to physical sector (that is to the subspace of constraints and gauge conditions) 
acquires the usual canonical form with the physical Hamiltonian contributed by the second 
term of (|5.12|) and H (q, p). Note that generically, especially for the so-called parametrized 
systems with H Q (q,p) = in canonical gauge conditions should explicitly depend 

on time, x M (<?) = X^{qi^)i in order to generate the dynamical evolution in reduced phase 
space theory [14], [J. Therefore, the reduced symplectic form generates a nontrivial contri- 
bution to the physical Hamiltonian, proportional to the time derivative of the embedding 
functions ( [5.6| ) explicitly depending on t, q = e(£, t). The total physical Hamiltonian then 
takes the form 



#red(£,7T) 



H (q,p) - p., 



dt 



(5.14) 



-I <?=<?(5.t).p=p(£:?0 



Canonical quantization of such a classical system runs as usual along the lines of a 
particular representation and operator realization in the Hilbert space of the theory. In 
the one-loop (linear in h) approximation with the Weyl ordering of the above Hamiltonian 
this quantization is basically exhausted by the unitary evolution kernel K(t, £|t', £') of the 
Schrodinger equation. In the coordinate representation it is given by the well-known 



Pauli-Van Vleck-Morette ansatz 15 



det 



2nh d£ A d£ B ' 



V 3 



(5.15) 



which involves the principal Hamilton function of physical variables S(t, £') - a clas- 
sical action evaluated at the classical extremal passing the points £' and £ respectively 
at initial t' and final t moments of time. The preexponential factor here is built of Van- 
Vleck determinant and guarantees in the approximation linear in fi the unitarity of the 
Schrodinger evolution of the physical states *ff(t, £) 



in the Hilbert space with a simple L 2 inner product 

(*l|* 2 )red = / d^KOMO- 



(5.16) 



(5.17) 



The unitary map between the reduced phase space quantization of the above type and 
the Dirac quantization of Sects. 2-4 consists in a special relation between the two-point 
kernel (|4.3|) with prefactor ( |4.9|) and the Schrodinger evolution operator ( |5.15| ) [p~3] , |^, |I1 . 
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This relation is based on the equality of the principal Hamilton functions in the original 
constrained theory and the reduced one and its corollary - the relation between their 
Van-Vleck matrices 



S(t,Z\t',e) = S(q,q': 



d 2 S(t,Z\t',?) 



Decomposing the gauge-fixed matrix (|4.10|) in the basis of vectors (e l A , V^) 
(|5.11|) one then easily finds the needed relation [JT3], |], [14 



K(t,£\t>,0 = const (j^f 2 K(q,q') (J^J 



1/2 



(5.18) 
(5.19) 
and using 

(5.20) 



q=e(t,t),q>=e(?,t>) 



This relation implies that the kernel K(q, q') similarly to the Schrodinger propagator 
K(t, £') can be regarded as a propagator of the Dirac wavefunction ^(q) in superspace. 
Indeed, introducing the following map between *&(q) and \I/(t,£) 



J\ 1 ' 2 



(5.21) 



and taking into account the relation (|5.7| ) between the integration measures on superspace 
and the physical space S, one finds that the propagation law ( |5.16|) in £-space can be 
regarded as a projection onto £(t) of the following propagation of the Dirac wavefunction 
^(q) from the initial Cauchy surface S(t') in g-space to the entire superspace 

¥(<?) = J dq' K(q, q') 5( X (q', t'))J{q', -dS(q, q')/dq') vfr(g') + 0(h). (5.22) 

Here the actual integration runs over the initial physical space S(t'). However, the in- 
tegration measure involves not just local quantities at this surface, but also the normal 
derivatives of the kernel (or the wavefunction *&(q') itself) arising in the one-loop approx- 
imation as a Hamilton- Jacobi argument —dS(q,q')/dq' of J(q,p). 



5.2. General gauge conditions 



Quantum reduction of the above type has been built in JIB , [TJ]] for a particular case 
of coordinate gauge conditions, x M (g, t) = 0, independent of momenta. The extension to 
general unitary gauges of the form ( |5.3| ) is less obvious, because such gauges no longer 
determine the coordinate physical space, but rather a (2n — m)-dimensional subspace in 
the entire pase space of (q l ,Pi). Therefore, the additional reduction of its dimensionality 
to 2{n — m) due to solving the constraints ( |1.2| ) leaves us with the physical phase space 
without a natural decomposition into coordinates and momenta. On the other hand, we 
need such a decomposition (choice of polarization) for sake of quantization in coordinate 
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representation. This decomposition can, however, be built if we model the p-dependent 
gauges by the coordinate ones and then repeate the construction of the above type. The 
idea is that semiclassically, the momenta equal the Hamilton- Jacobi values up to quantum 
fluctuations, p = dS/dq + 0(h), so that the coordinate gauges can be obtained from the 
general unitary gauges by the replacement 

OS 

X*(q,P,t) -> x^q^t) = x M (<7, -Qj,t). (5.23) 

Then the reduction procedure is rather straightforward, although it involves an addi- 
tional step of the unitary transformation from the physical variables in the initial gauge 
X^{q,p) = to those of its coordinate version ( |5.23| ). 

The peculiarity of the reduction in p-dependent gauges is that, in contrast with ( |5.12| ), 
the canonical transformation of the original symplectic form Pidq 1 to the physical one 
7r A d^ A is not contact. It involves a nontrivial generating function F(q,£,t) 

(pidq'-Hodt) u ^ = Ti A di A - H rcd dt + dF(q,£,t), (5.24) 

so that the original action fll.lj) in terms of physical variables acquires extra surface terms 
at the final and initial moments of time t±| 

S[q,p] MT = f t+ dt(7c A £ A -H tcd ) + F(q + ,Z + ,t + )-F(q_,Z_,t_). (5.25) 

The actual form of the generating function F(q,£,t) depends on the choice of gauge 
conditions and the choice of physical variables (£ , tt^). It can be rather complicated, 
but, fortunately, we do not need it explicitly. 

Important property of the surface terms in (|5.25 ) is that they take care of the bound- 



ary conditions on the histories (q(t),p(t)) at t±. Natural boundary conditions for the 
variational problem with the action ( |1.1|) - left-hand side of eq. (|5.25| ) - correspond to 
fixed coordinates 

q\t ± )=qi (5.26) 

and arbitrarily varied momenta. On the contrary, the variational problem for the reduced 
phase-space action - the integral term in the right-hand side of fl5.25Q - assumes fixed 
physical coordinates 

i A {t ± ) = (5.27) 

but the reduction from (q,p) to (£, 7r) intertwines coordinates with momenta, q = g(£, tt). 
Generating functions in the surface terms of the right-hand side of ( |5.25| ) reconcile these 
two different boundary value problems. 

2 For notational reasons in this section we supply the arguments of two-point quantities by ± instead 
of primed and unprimed labels above. 
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Another important aspect of the boundary conditions (|5.26|) is that the values q± can- 
not be arbitrary. In complete analogy with eq.( 5.18| ) they should belong to the physical 
subspace of the full phase space, that is to satisfy both the constraints and gauge condi- 
tions. This property was obvious in coordinate gauges which directly restricted possible 
values of q± to the surface ( |5.6| ) of gauge conditions. For p-dependent gauges one could 
think that the flexibility of choosing the momenta at t± completely releases the values q± 
to be arbitrary. This is, however, not the case. One way to see this is to calculate the 
action at the extremal satisfying classical equations of motion and boundary conditions 
(|5.26|) . The resuling principal Hamilton function S(q + ,q-) uniquely determines the end- 
point momenta, p± = ±dS /dq±, and the end-point gauge conditions take the form of the 
coordinate gauges of the form ( |5.23|) 

XUQ±) = X»(q±,±dS/dq ± ,t ± ) = 0, (5.28) 

that determine the embedding of the following two subspaces into the coordinate space 

E± : qi = e l ( V i,t ± ), ^(e'fo, t±)) = 0. (5.29) 

Here the internal coordinates r\ A no longer coincide with the physical coordinates £ A of 
the reduced action. 

Thus, the principal Hamilton function reduced to subspaces (|5.29| ) takes the form 
5 (<?+;<?-)| 9±=e(j7±)t±) = S(t+,r}+\t-,Tl-), (5.30) 

«s(t + ,7 7+ |t_,7 7 _) = ^ + ,e + K-,e-)+^+,e+,i + )-^- ) e-^-), (5.31) 

where S(t + , £_) is the principle Hamilton function of the reduced system in the 

p-dependent gauge with fixed end-point coordinates ( |5.27| ) and 

F(v±,t±,t±) = F(e( V± ,t±),£±,t ± ) (5.32) 

represent the generating functions of canonical transformations from £± to i]± at t±. The 
whole expression on the right-hand side is a function of rj±, the coordinates £± being 
expressed in terms of T)± as solution of equations 

dS(t +t £ + \t-,Z-) ^ dF( V± ,Z ± ,t ± ) _ n 

(these in turn follow from the corollary of eq. fl5.24 ), tta = —dF(q,^,t)/d^ A , and the 
Hamilton- Jacobi value of the physical momentum, ir± = ±dS(t+,£+\t-,£-)/d£±). 

In this gauge we can now repeate the construction of the Dirac two-point kernel of 
Sect. 4. We identify the covectors (XiiXv) with gradients of coordinate gauge condi- 
tions, dx±{q)/dq l . Then we use fl5.30p and again have the relation ( |5.20 ) with the kernel 
K (£,£!*',£') replaced by 



£(*+>»7+l*-,»7-) 



' i d 2 S(t + ,r] + \t^,r]_) 
2irh drj + drj_ 



1/2 i_ 



*S(t + , T)+\t-, ri -) 

> (5.34) 
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and the factors (J, J') = J± (|4.11| ) calculated in the coordinate gauges of the above 
type (f>. 28|) . Apriori, they differ from the Faddeev- Popov determinants corresponding to 
p-dependent gauges, but as shown in Appendix A, eq. ( |A.7|) , 



x' 



{x^q,p)Mq,p)} , (5.35) 

p = dS j dq 

and all the J/f-factors coincide with true Faddeev-Popov matrices in generic unitary 
gauges. 

The interpretation of the kernel ( |5.34|) is obvious. It is built of the principal Hamilton 
function ( |5.31|) corresponding to fixed boundary values of //-variables at t±. As dis- 
cussed above, it contains as total derivative terms the generating function F(r?,£,t) of 
the canonical transformation from £± to r]±. At the quantum level, in the semiclassical 
approximation linear in h, the kernel of the corresponding unitary transformation reads 



U{t,v\t) 



det 



2nh drj dt; 



(5.36) 



If we represent the kernels of the above type as matrix elements of the corresponding 
unitary operators in the coordinate (£ and 77) representations 



K(t + ,Ti + \t-,TjJ) =< rj + I fc(t+,tJ) 1 77- >, 

u(t, v \o = <r] \u(t)\z>, 

K{t + , C+\t-, e_) =< e+ I K(U, t.) \C->, 



(5.37) 
(5.38) 
(5.39) 



then the kernel in question ( |5.34j) reads as follows. It is nothing but the unitary evolution 
operator of the reduced theory in generic p-dependent gauge K(t+,t-) (whoes kernel is 
defined by eq.(5.10) with the principal Hamilton function subject to fixed (£,£') = £±) 
unitarily transformed to the representation of 77-variables 



(5.40) 



This unitary equivalence relation can be directly checked by the composition of kernels 
(|5.36| ) and ( |5.15| ) and calculating the corresponding inegrals over £± by stationary phase 
technique^. 

This reveals the role of extra unitary transformation arising in the reduction to physical 
sector in generic momentum-dependent unitary gauges ( |5.5| ). A similar mechanism of 
modelling the non-unitary (relativistic) gauges, which involve the Lagrange multiplyers 
and their derivatives, by unitary ones was recently proposed in where the Dirac 
two-point kernel was obtained by a direct calculation of the one-loop path integral in 
gauge theory. The construction of such gauges is based on the knowledge of solutions of 
linearized wave equations in a corresponding relativistic gauge. 



The contribution of stationary points - solutions of eqs.(5.33) - yields the prexponential factor of 



(5.34), while its phase is achieved as a linear combination of Hamilton- Jacobi phases (5.31) 
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6. The physical inner product: Hermit icity and gauge 
independence 



The auxiliary inner product ( |2.8| ) cannot serve as an inner product for physical states 
because it is not well defined. In view of quantum constraints the physical states have 
a distributional nature | = "<5(T)" | ^ aux ) with somehow determined m-dimensional 
delta-function of non-abelian operators T M , and their naive bilinear combinations are 
divergent because ["<5(T)"] 2 ~ <5(0)"<5(T)". At most the auxiliary vectors | V&aux) partici- 
pating in the construction of the physical states can be required to be square-integrable 
in L 2 sense and thus induce a finite inner product for | ty) (which is the idea of the so- 
called refined algebraic quantization of constrained systems fT8| , pjj). Another approach 
consists in the unitary map from the Dirac to the reduced phase space quantization of 
the previous section. Reduced theory has a trivial inner product ( |5.17|) which induces a 
correct physical product in the Dirac quantization scheme 

(* 2 |*l) = (*2^l)red (6.1) 



on account of the relation ( |5.21| ) and the change of integration variables ( |5.9| ). The result 



for semiclassical states of the form ( |4.1| ) looks like |2(], [13], [5|] 

(¥ 2 |*i) = / dq^l{q)8{ X {q))J{^dS/dq)^ 1 {q)+0{h). (6.2) 
Here x(<?) = X M (<?) is a set of gauge conditions delta-function of which 

S(x) = I[Wq)), (6.3) 

determines the (n — m)-dimensional physical subspace S embedded in superspace and 
J(q,dS /dq) is a corresponding Faddeev- Popov determinant^. 

This explains the nature of the physical inner product in the Dirac quantization. Its 
measure contains the Faddeev-Popov determinant which depends on the choice of the 
constraint basis and transforms under the transition ( |1.7|) to another basis as a density 
of weight 1 in the space of gauge indices (^4.12|) as compared to -1/2 weight of the Dirac 



wavefunctions (|3.8j ). This proves the invariance of the physical inner product under this 
transformation and shows that at the quantum level it is not only canonical but also 
unitary. Our purpose now, till the end of this section, will be to discuss the Hermiticity 
properties of physical observables relative to this inner product and the gauge indepen- 
dence of their matrix elements. 



4 In eq. (|6.S|) the product S&^i = ^2* -Pi exp[i(Si — Sy/Ti] involves two different Hamilton- Jacobi 
functions, so that it seems ambiguous on which Lagrangian manifold (p = dSi/dq or p — dS2/ 'dq) the 
relevant momentum argument of J(q,p) should be constructed. One should remember, however, that in 
semiclassical expansion the integral is calculated by the stationary phase method in which a dominant 
contribution comes from the stationary point satisfying dSi/dq — dS2/dq. This makes these Lagrangian 
surfaces to coincide in the leading order, their difference being treated perturbatively as expansion in H. 
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The semiclassical physical inner product (|6.2|) can be rewritten as an auxiliary inner 
product of physical states with a nontrivial operatorial measure 

(*i|* a ) = J8(x) I * 2 > + 0(h). (6.4) 

Here the operator ordering in operators of the ghost determinant and gauge conditions 
is unimportant because it effects the multiloop orders 0(h) that go beyond the scope of 
this paper. The matrix element of the physical observable Oj is therefore 

(*i| 6 l |* 2 ) = (*a| 6j* 2 ) = ( *x| J5(x)Oi\ * 2 ). (6.5) 

To check the Hermiticity of Oi we have to show that the expression 

( 6 I ^ 1 \ tf 2 ) - a,* 2 ) = ( 0\J5(x)\ * 2 ) - ( *i| JS(x)dj\ * 2 ), (6.6) 

where a dagger denotes Hermitian conjugation with respect to the auxiliary inner product, 
is vanishing. From ( |2.13| ) it follows that 0\ = Oj—ihUj X +0(h 2 ) (we consider the algebras 
of observables with Ufj = 0), and the expression above takes the form 

( *x| [ Oj, J5(x) ] - ihU} x | * 2 ) = 0(h 2 ), (6.7) 

which, as shown in Appendix A, is vanishing in the one-loop approximation. Thus, 
the physical observables are semiclassically Hermitian with respect to the physical inner 
product (|6.2|). 

Another important property of this product and matrix elements of observables is their 
independence of the choice of gauge conditions participating in their construction. 

The gauge independence of the inner product itself is based, as shown in |5|, ^T|], on the 
fact that it can be rewritten as an integral over (n — m)-dimensional surface £ of certain 
(n — m)-form which is closed in virtue of the Dirac constraints on physical states 

(* 2 |*i) = X^^^' dujin ' m) = °- ( 6 - 8 ) 

It follows then from the Stokes theorem that this integral is independent of the choice of 
S or equivalently of the choice of gauge conditions specifying the physical subspace. The 
form uj( ri ~ m '> in the one-loop approximation equals 

^ n ~ m) = rfgn ^l A ^f"^ n...,^ 2 Vi— +1 -..V^ *x, (6.9) 

and its closure is a corollary of the continuity equation ( |4.5|) for ^2^1 

^-(v;*^) = t^*;*!. (6.io) 

The evaluation of the matrix element of the observable O semiclassically involves the 
evaluation of the quantity 

(* 2 | 0,*!) = / cu^O^q, dSi/dq) + 0(h), (6.11) 
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which will be also gauge independent provided the continuity equation holds for the 
quantity ^l(q)^fi(q) Oj(q,dSi/dq). But this equation will also be a corollary of ( p.lO| ) 
because, as shown in Appendix A, on the Lagrangian manifold of phase space 



^dq* 



dq 



P 



= (6.12) 

dS/dq 

in view of gauge invariance ( [2.10| ) of observablesf]. Thus, as it should have been expected 
from the theory of gauge fields [17|, |16[ the gauge independence of the physical matrix 
elements or expectation values of observables follows from the gauge invariance of the 
latter. This is true not only at the formal path-integral quantization level, but also in the 
operatorial Dirac quantization scheme. 



7. Conclusions 

Thus we see that, despite a complicated non-abelian nature of the formalism, Dirac 
quantization of generic constrained systems is remarkably consistent and reveals rich ge- 
ometrical structures beyond the lowest order semiclassical approximation. Geometrical 
covariance of operators and physical states takes place not only in the coordinate config- 
uration space of the theory, but also in the space of gauge transformations. The quantum 
formalism turns out to be covariant in the sense of unitary equivalence with respect to 
generic symmetries of a classical theory including arbitrary change of the constraint basis. 
The operators of physical observables turn out to be Hermitian with respect to a physical 
inner product and their matrix elements are gauge independent in full correspondence 
with similar properties in the classical domain. All these properties were obtained pertur- 
batively in the one-loop approximation of semiclassical expansion but, no doubt, they can 
be extended to multi-loop orders, though, apparently by the price of growing technical 
complexity. 

It should be emphasized that this remarkably general picture of quantum invariances 
was obtained in the Dirac quantization of constrained systems, when all quantum con- 
straints are imposed on physical states. These states can be regarded a quantum trunca- 
tion of the BFV(BRST) quantization with CP-ordered form of the nilpotent BRS opera- 
tor in the extended relativistic phase space of original (q, p) and ghost canonical variables 
(C, V) 0, H OH- As mentioned in Introduction, the main motivation for such a setting is 
the quantization of zero modes of extended objects for which another (unitary inequiva- 
lent) quantization in the Fock space based on the normal Wick ordering in the ghost sector 
is not applicable. The latter is usually applied to oscillatory modes of field theoretical 



5 In eqs.( |6.9[) -( |6.1l|) the difference between the momentum arguments reiated to different Hamilton- 
Jacobi functions p — dSi/dq and p = dS^jdq should also be treated perturbatively in h and, thus, goes 



beyond the one-loop approximation (see the footnote to eqs.(S.2)-(|0|)) 
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models, including strings and low- dimensional CFT This quantization after trunca- 
tion to (g,p)-sector results only in a half of the initial constraints imposed on physical 
states (and in view of this admitting the central extension of quantum algebra). For 
them the geometric properties of the above type are much less known if at all available at 
such a general level not resorting to a concrete harmonic oscillator decomposition of fields. 
Thus it seems interesting to try extending the above geometric methods to this quantiza- 
tion scheme, although it is very hard to expect that the general coordinate invariance in 
curved configuration space of the theory can be fruitfully combined with the Fock space 
representation. (The latter, as is well known, is basically applicable to linearized fields 
belonging to the tangent vector space rather than to the configuration space manifold.) 

Another difficulty with the field theoretical (infinite dimensional) extension of the 
above technique consists in ultraviolet infinities and the accompanying them anomalies 
that would generally violate gauge symmetries at the quantum level. As is well known, 
this problem cannot be analyzed at the level of generality adopted in this paper, that is 
for gauge theories of general type. A usual approach would consist in applying the above 
formalism to a given field model within a particular regularization scheme that would 
render all the quantities finite, singling out the potential anomalies of gauge symmetries 
and using the results as selection criteria for viable systems. Like in string models this 
might lead to restrictions on the number of fields (dimensionality of target spacetime) or 
other parameters of quantized model. More detailed discussion of this issue can be found 
in 



As far as it concerns the present results, their direct implications can be expected in 
quantization of zero modes in various field theoretical models or extended objects. In 
quantum cosmological context the Dirac quantization scheme is implemented in the sys- 
tem of Wheeler-DeWitt equations. The minisuperspace (zero) mode of the full superspace 
of quantum cosmology cannot be quantized in the Fock representation, because it has a 
ghost nature p5, |26|, 27], so only the coordinate representation and relevant methods 
of this paper can be used for the studies of the quantum Cauchy problem for the early 



inflationary Universe |14|, |23|. The gauge invariant observables in early cosmology, the 
theory of which is very important for understanding the formation of structure |28| and 
back reaction phenomena, should be analyzed by the technique of the above type. In this 
respect, especially important becomes the reduction technique in generic unitary gauges 
that should be applied in the theory of cosmological perturbations of |2Bj in order to 
study the quantum back reaction effects in effective equations of inflationary dynamics 
|[29|| . The potential range of applications is even wider in D-brane dynamics of string 
theory |Kj and brane cosmology ][H|] . Zero modes of D and p-branes in nonperturbative 
string theory or brane- worlds in cosmology jj^] , separation of their dynamical from purely 
gauge properties, description of brane dynamics in the bulk in the formalism of quantum 
Dirac constraints and many other issues are subject to the methods of this paper. 
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A. Hermiticity of observables 



To prove eq.( |6.7| ) note that 

( * x | [ 6j, J5(x)}\ * 2 > = ihjdq{ Oj, J5( x )} *t* 2 + 0(h 2 ). (A.i) 

The Poisson bracket commutator here can be transformed by using the cyclic Jacobi 
identity and the weak gauge invariance of constraints ( |2.10| ) 



dq{0 I ,J6( X )}n*2= I dq5( X )J [ + J~Y{T^ { X v , Oj}} \ *l* 2 

dx» 

The last term here can be integrated in a special coordinate system on superspace defined 
by eqs. (|5.5| ) and 



(A.2) 



J dq ^ J {Oj, X*} *t*2 = ~jdq Six) M 



dx 



d 

de» 



(A.3) 



To transform this expression further let us derive several useful identities. First of all, 
note that any function of phase space variables f(q,p) when restricted to the Lagrangian 
manifold defined by the Hamilton- Jacobi function S becomes a function on superspace 
f(q,dS/dq). Its derivative with respect to coordinate M (of the new coordinate system) 





ovi' v < )( i , 

can be simplified to 



d 



e , (df_ + 9f d 2 s • 

M \ dq { dpk dq k dq 1 



(A.4) 



(A.5) 



in view of eq. (|5.10| ) and the differentiated version of the Hamilton- Jacobi form of con- 
straints 



d 2 s 



dT„, 



(A.6) 



M dq 1 dq k dq k 

Similar derivation shows that the gauge transformation of this function on the Lagrangian 
manifold, generated by the vector flow V* , coincides with the Poisson bracket of f(q,p) 
with the constraint 



_d_ 

' r dq i 



f k 



dtT 

dq, 



(A.7) 



evaluated certainly at p = dS / dq. With these identities and using the derivatives of the 
measure M 



dM 
~d¥ 



del 

-M — - 



dM 



dq i ' d^ A 



-M 



dq 1 



(A.8) 
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(e* and e l A are defined by (|5.S| )), the relation ( |5.10|) and the Jacobi identity for Poisson 



brackets one can obtain the following gauge derivative 



d 

m 1 



j_ 

M 



J_j-iu(dVl 
M " I dq 1 




(A.9) 



Using this equation for J/M and applying (|A.6|) and (|A.7|) to other quantities in ( |A.3|) 
we finally obtain 



dqS( X )M 



d 



+ 



+ 



J dq 5(x) J.rY{T„ { X u , Oj}} *t* 2 

" d 



dq 



M 9 



J_ 



(A.IO) 



The second term here is vanishing in view of the continuity equation (or more precisely 
0(h) for different SE^ and Sf^)- The third total derivative term is vanishing in view of zero 
boundary conditions for physical wavefunctions at the infinity of the physical space E. 
Collecting equations (|A.2|) , ( |A.3j ) and ( |A.10| ) together we see that the matrix element of 



the commutator reduces to the one term containing the trace of the structure function Uj 



which exactly cancels out in the equation (577). This proves the Hermiticity of observables 
in the physical inner product. 



B. Abelianization procedure on Lagrangian manifolds 

It is well known that when the structure functions of the Poisson-bracket algebra 
of constraints (|1.3j ) are not constants this algebra is open ||: the commutator of two 



consequitive transformations S^f = {/, T M } of any function on phase space 

{{/, T M }, T v } - {{/, T v }, T M } = {/, T\} + {/, U^} T\ (B.l) 

is a linear combination of these transformations only on the constraint surface T\(q, p) = 0. 
Semiclassically, the restriction to this surface takes place on the Lagrangian manifold of 
phase space 

» - 1 < B - 2 > 

defined by the Hamilton- Jacobi function of a semiclassical state, satisfying the Hamilton- 
Jacobi equation ( f4.4|) . As was shown in Appendix A, the action of the constraint gener- 
ators on f(q,p) (in the sense of Poisson brackets) on this surface ( |A.7|) can be generated 



by directional derivatives along a special set of vectors (|4.6|) . Therefore these vectors 
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can be regarded as gauge generators on the Lagrangian manifold of phase space. They 
have a closed Lie-bracket algebra 

- Vt-^ = £/lV$ (B.3) 
and, therefore, can be abelianized p4| by recombining these vectors with the aid of the 



matrix of the inverse Faddeev- Popov operator constructed out of some admissible gauge 
conditions = V*dx M /<9<f 



0. (B.4) 



Abelian generators alow one to construct preferred parametrization of the coordinate 
manifold ( |5.5| ) with a special internal coordinates on a physical space E. Note that the 
new coordinates £ A in ( |5.5| ) as functions of the original coordinates are not necessarily 
gauge invariant. Abelianization procedure can render them invariant as follows. Demand 
that the reparametrization functions in ( |5.5|) satisfy the equations 

which are integrable in view of the abelian nature of TZ l . Then the identity e l d£ /dq l = 
implies the gauge invariance of £ (q) 

dq i 

Many equations above simplify with this preferred parametrization of physical space E. 
Indeed, eq.( [B.5|) implies that = in the equation (|5.10|) and the equation (|A.9|) for 
the factor J/M performing the map from the Dirac quantization to the reduced phase 
space quantization takes the form of the continutity equation for a one-loop prefactor 
d_ 
dq l 

which means that this factor with the measure M constructed in this parametrization 
represents a kinematical solution of the continuity equation. 



— (V* J/M )=U^J/M, (B.7) 
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